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We study the cyhnder to sphere morphological transition of diblock copolymers in aqueous so- 
lution with a hydrophobic block and a charged block. We find a metastable undulated cylinder 
configuration for a range of charge and salt concentrations which, nevertheless, occurs above the 
threshold where spheres are thermodynamically favorable. By modeling the shape of the cylinder 
ends, we find that the free energy barrier for the transition from cylinders to spheres is quite large 
and that this barrier falls significantly in the limit of high polymer charge and low solution salinity. 
This suggests that observed undulated cylinder phases are kinetically trapped structures. 
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Diblock copolymers, formed by joining two immisci- 
ble polymers at one end, are a model system to study 
many aspects of self-assembly and order-order transitions 
from one self- assembled state to another 1, 2, 3j. Self- 
assembly in diblock copolymers is driven by two com- 
peting effects: the immiscibility of the two blocks drives 
their phase separation, yet the connectivity of the blocks 
prevents a true macroscopic phase separation. Even in 
the simplest case of copolymer melts, this competition 
leads to the formation of a variety of generic structures, 
including spherical and cylindrical micelles, lamellae and 
bicontinuous gyroid phase 0, |^ . The selection of these 
various morphologies is determined by a balance of sur- 
face tension between microphase separated regions and 
the chain stretching of the polymers and can, in certain 
cases, be analyzed purely from the point of view of the 
interfacial geometry 6]. 

In selective solvents, on the other hand, the morphol- 
ogy of copolymer structures is also affected by solvent- 
polymer interactions, widening the both the range of pos- 
sible structures and the ability to control morphological 
transitions by tuning the solvent conditions. Moreover, 
in solution, the volume-filling constraints of the molten 
state are relaxed. In this way, diblock copolymers are 
the analogues of surfactants, and exhibit many of the 
typical self-assembled morphologies: spherical micelles, 
wormlike (cylindrical) micelles, and lamellae 0. Despite 
the similarities, however, these systems also display a va- 
riety of novel structures, including wormlike micelle net- 
works 1^ [lO| , large compound vesicles [HI, •> aiid 
undulated cylinders 0, 0, [TO, [ll . 

In this paper, we study the undulated cylindrical mor- 
phology by considering axisymmetric, diblock copolymer 
micelles dispersed in aqueous solution, concentrating on 
the case of copolymers composed of a hydrophobic block 
and a polyelectrolyte block. In part, this focus is due to 
the fact the self-assembly in polyelectrolyte copolymer so- 
lutions has been well characterized experimentally |T^ . 
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While similar experiments have address the morphologi- 
cal transition between wormlike micelles and spheres for 
neutral block copolymers |0, ^3 5 added control pro- 
vided by modulating the charge and salt conditions of 
the polyelectrolyte domains allow for a systematic in- 
vestigation of, in particular, the cylinder to sphere (CS) 
transition. Bendejacq et al j3| have performed such ex- 
periments at controlled salinity and ionic strength and 
find that when the polyelectrolyte domains are highly 
charged, the CS transition proceeds by the undulation of 
the cylindrical micelles followed by the eventual breakup 
into spheres. At certain intermediate ionic strengths and 
salt concentrations, however, undulated cylinders were 
also observed to have very long lifetimes, and were sug- 
gested as a candidate for a new intermediate phase be- 
tween cylinders and spheres in PS-PA A diblock copoly- 
mers. 

We develop a mean-field model for diblock copolymer 
micelles with interfacial geometry between straight cylin- 
ders and individual spheres. Experimentally, this can 
be realized by first forming the cylinders, then changing 
solvent conditions to frustrate the geometry toward the 
formation of spherical micelles, as was the method in the 
experiments of ref. 13]. To this end, we construct a class 
of micellar morphologies whose geometry can be contin- 
uously varied from cylindrical to spherical and compute 
the free energy within a variational, mean-field scheme. 
Our main result is that under certain salt and charge con- 
ditions, cylindrical micelles become unstable to the for- 
mation of undulated cylinders. These undulated cylin- 
ders are metasta ble. however, and their observation by 
Bendejacq et al. [T^ must indicate a large kinetic bar- 
rier to the breakup of undulated cylinders to spheres. 
Our model allows us to calculate an upper bound on the 
free-energy barrier for fission of a cylinder into spheres, 
and we find a strong dependence of the barrier height on 
salt concentration of the solution and the effective charge 
of the polyelectrolyte blocks. Combined these calcula- 
tions enable us qualitatively to understand the "phase 
diagram" observed in reference 13]. We also discuss re- 
sults on the shape of the micelle ends. 

Our approach is to construct the simplest possible 
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model for the energetics of undulated cylinders that con- 
tains the essential physics. In this spirit, we will ig- 
nore the fluctuations of the chains around their prim- 
itive paths and assume the chains in the corona are 
equally stretched. Strictly speaking, this analysis pro- 
vides only an upper-bound on the true free energy. Nev- 
ertheless, similar calculations are known to be very useful 
in understanding the equilibrium morphologies of diblock 
copolymer melts 0. Moreover, for charged polymers 
in particular, which tend to be strongly-stretched, these 
techniques have been developed extensively by Borisov, 
Zhulina and collaborators 0. Recently, Jusufi et 
al. 18] have performed extensive simulations and detailed 
analytical analysis on charged star polymers and found 
excellent agreement between simulations and a mean- 
field electrostatic theory. 

The paper is organized as follows: In section ID we dis- 
cuss the details of our model for the diblock copolymers. 
In section^ we present phase diagrams and results on 
a metastable undulated cylinder phase. We also discuss 
the shape of the cylinder ends and compute estimates for 
the free energy barriers. In section IIIIl we discuss our 
results. 
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FIG. 1: The unduloid family of surfaces for some particular 
values of e, spanning shapes from cylinders to a necklace of 
spheres. These surfaces have constant mean curvature given 
by (2i?c)~^, where Rc is the radius of the cylindrical state at 
e = 0. 



I. MODEL FOR UNDULATED CYLINDERS 

As a model for the interface we consider the unduloid, a 
family of constant mean curvature surfaces that appear 
in the study of soap films These have also been 

used to describe the pearling instability in fluid mem- 
brane tubes 20]. For our purposes, the unduloid surfaces 
have the convenient property that their geometry can be 
tuned continuously and rather simply from cylindrical to 
spherical. This family of surfaces is described by two pa- 
rameters: e which varies from for the cylinder to 1 for 
a necklace of spheres, and the size scale Rc which sets 
the radius of the cylinder at e = 1 Intermediate 
values of e describe an undulated cylinder surface (see 
Fig. nj. Following ref. |20|, the surface can be repre- 
sented parametrically with a vector that maps li, the co- 
ordinate along the cylinder axis, and (/), the azimuthal co- 
ordinate, to points on the interface between the charged 
and hydrophobic blocks, X('U, (j)) = z{u)z + r{u)r, where 

^ . 1 + e sin t , , 

= Rc dt . (1) 

J-7T/2 V 1 + 2e sm t-^ 

and 

r{u) = + 2esinii + e2. (2) 

We see that X(ii,(/)) is periodic diS u ^ u -\- 27r. It is 
important to note that, while Rc describes the radius of 
the cylinders, the radii of the spheres at e = 1 is 2Rc. 
Note that at small e, the unduloid describes a simple 
sinusoidal modulation of the cylinder radius, 

r{z)/Rc = 1 + e cos{l /Rc) + O(e^) . (3) 

In the case of soap films, unduloids are necessarily un- 
stable due to their finite length. For diblock copolymers, 
however, we will see the situation is far more complicated. 

The free energy of polyelectroyte copolymer micelles is 
a sum of four terms, 

F = Fc^st + Fint + Fb^st + Fci, (4) 

where Fc^st is the stretching energy of the chains in the 
core, Fint is the interfacial energy between the core and 
dissolved brush domains, F^j^gt is the stretching energy of 
the brush, and Fd is the contribution to the free energy 
from the charged groups on the chain and the ions in 
solution. In the remainder of this section we provide 
results and motivation for the expressions, leaving their 
detailed evaluation to the appendix. 

A. Micelle Core and Interface 

If the polymer block composing the core region is suffi- 
ciently hydrophobic, then the concentration of solvent in 
the core region is negligible and the chains in the core are 
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FIG. 2: On the left, a schematic of the structure of the strongly- stretched corona. The corona chains are approximated by the 
equal stretching hypothesis of Gaussian chains along the average trajectories (dashed). On the right, the chains in the corona 
are normal to the interface surface for f3 = 1 and radial for = 0. The density of chains is indicated by the density of primitive 
paths. The chain density is least in the undulated necks because fewest chains reach the interface there from the core region. 



in the molten state with constant density po. Thus, the 
techniques of self-consistent strong-segregation theory of 
diblock copolymers in a melt can be applied 5]. In the 
strong- segregation limit, the width interface between do- 
mains is small compared to the size of the domains, and 
the free energy can be expressed in terms of the geome- 
try of the interface. To this end we employ the "wedge 
model" of Milner and Olmsted in which we partition 
the core region-contained by the unduloid surface-into a 
series of wedges small enough so that all chains within 
effectively experience a single, fixed geometry. The total 
stretching free energy is computed by, 

where Jy dV denotes an integral over the volume of the 
core domain, appropriately apportioned into infinitesimal 
wedges, he is the distance along the chain trajectory to 
the interface, po is the density of monomers in the core, 
Nc is the molecular weight of the chains in the core and 
a is the Kuhn length. We assume that the chains meet 
the interface along the surface normal. One advantage of 
using the unduloid representation of the interface rather 
than a cylinder with a cosine modulation of the radius 
is that, for all values of e, the chain paths do not cross, 
allowing us to follow the transition all the way from the 
cylinder morphology at e = to a chain of spheres at 
e = 1. 

We show in the appendix that the stretching energy 
per chain can be put into the form 

F.^st _ ttRI S{e) 
nksT mco? V(e) ' ^ ^ 



where n is the number of chains, V(e) = V/ii'^R^ is the 
dimensionless volume of the micelle core, V is the core 
volume, and 5(e) is the dimensionless stretching energy 
for the chains in the core. Since this term represents 
the entropic cost of confining a chain within a brush, the 
free energy simply scales as R^/NcO^^ where NcO? is the 
unperturbed radius of gyration of the polymer block. 

At the interface between the core and dissolved brush 
domains, the hydrophobic block is exposed to solvent, 
as well as the polyelectrolyte block. Assuming that the 
interface is sufficiently narrow, we can capture energetic 
cost of this exposure by introducing a free energy per unit 
area, or surface tension, 7. In general, such a surface ten- 
sion will have a complicated dependence on the relative 
concentrations of solvent and polymer at the interface, 
free charge and salt concentration, etc. We will treat, 7, 
as an unknown phenomenological parameter effectively 
fixed for a given system at constant temperature. The 
interfacial energy per chain is then given by 

F^nt ^ lA{e,R,) ^ 2-iN,A{e) 
nksT n poRcV{e) ' ^ ^ 

where 7 is the surface tension (in units of ksT)^ A(e, Rc) 
is the interface area, and A{e) = A{e^ Rc)/7t'^ R^ the di- 
mensionless area of the unduloid surface. 



B. Micelle Corona 

In the limit that the concentration of free ions (both 
counter ions and salt ions) in the brush is large, the pri- 
mary effect of these free charges is to maintain electroneu- 
trality 16]. In this limit, the so-called "osmotic-brush" 



4 



regime, the stretching of the chains is balanced by the en- 
tropy of the counterions confined to the brush, and the 
counterions locally neutralize the charge on the chains 
|2l| . In the present calculation, we consider chains that 
are equally-stretched 22], and neglect polymer trajec- 
tory fluctuations. The latter approximation can be jus- 
tified in the low-salt limit and the limit of large poly- 
electrolyte block monomer number, A/'^. Fluctuations 
of chain path around the mean trajectory have charac- 
teristic size, N^^'^a, while simple scaling considerations 
show that in the low-salt limit the brush height grows 
as Hb ~ N^a Thus, in the asymptotic, large 

limit, such fluctuations can be ignored. Further, Jusufi 
et al. have found that these these approximations 
lead to reasonable quantitative agreement with simula- 
tions for star polymers because the osmotic pressure of 
the counterions lead to highly stretched chains. How- 
ever, their model includes two effects which we neglect 
here: (1) the condensation of counterions, which we sim- 
ply absorb into a redefinition of the fraction of charged 
monomers, a, and (2) the impartial neutralization of the 
brush due to the escape of counterions to solution. 

A schematic of the chain model is shown is figure El 
We assume a brush height, and Gaussian chains with 
stretching energy 



hO?' J 



dA Hr 



where for simplicity we take the Kuhn length of the poly- 
electrolyte block to be, a, as it is for the insoluble block. 

To allow for appropriate generality, we introduce into 
our ansatz another parameter /3 which describes the av- 
erage direction of the chains in the brush. 



/3n + (1 - p)r 



(9) 



v//32 + (1-/3)2 + 2/3(1 -/5)n-f' 

where h is the unit normal of the interface (see Fig. [2j. 
In the regions of high negative mean curvature, the den- 
sity of chains, and hence the osmotic pressure, is much 
larger than the rest of the brush. Since charged brushes 
do not typically overlap [2^, we can allow the configu- 
rations to relax away from overlap in two ways. One on 
hand, the chains trajectories can tip radially away from 
the high density region, as is described by a /3 = 1 con- 
figuration. We could further allow the brush to compress 
by imposing a sinusoidal variation of the brush height, 
Hi){u) = Hi){l^vsm.u). For configurations without over- 
lap, we set = 0, and otherwise, we choose the smallest 
V such that the brush does not overlap. Indeed, we find 
over a large range of parameters the free energy favors 
uncompressed states, with v = {). By suitably rescal- 
ing Fh^st^ the dimensionless corona stretching per chain, 
B[Hb{u), e], can be defined as 

F.^st 3Rl B[n^,e,(3] 



nki 



(10) 



where Hb = H^/Rc is the height of the brush measured 
in units of the core size. 



Finally, we assume that each chain has a fraction 
a of charged monomers equally distributed along the 
chain. That is, we assume that the charge distribution is 
quenched. The counterion contribution to the free energy 
is then [l^ 



c+(ln(c+po)-l) (11) 
+c_ ( In (c_po) - l) - ( \n{csPo) - l) 



where c± is the position-dependent concentration of pos- 
itive and negative ions, Cg is the concentration of ions in 
bulk solution, and the integral is over the volume of the 
brush. The mean-field counterion densities should, prop- 
erly, come from solving the Poisson-Boltzmann equation 
for the charge distribution described by ac^(r), where 
Cm(r) is the local concentration of monomer in the brush 
which depends on the brush geometry (see Appendix). 
To derive this expression, we have assumed local elec- 
troneutrality, c+ — c— = acm which is valid in the 
osmotic regime . From the Poisson-Boltzmann equa- 
tion, we also have the relation c+c_ = for a 1 : 1 
electrolyte. This allows us to derive expressions for the 
ion concentrations in terms of Cm and Cg: 

acmir) 1 



c±(r) = + i VK-W)' + (12) 

Strictly speaking, we are considering the limit of the 
Poisson-Boltmann theory in the limit that the Debye 
screening length, is small compared to brush di- 

mensions 21]. This is valid for the case when the con- 
centration of free charge (counterions and salt) in the 
brush is sufficiently large. We note that the expression 
for the osmotic pressure used in a similar mean-field cal- 
culation for polyelectrolyte copolymers by Borisov and 
Zhulina [2^ can be obtained by taking the Cs ^ Cm limit 
of the above expressions. 

In general, one expects that some fraction of the coun- 
terions will be condensed onto the chains - we can ab- 
sorb this into a redefinition of the fraction of charged 
monomers a and assume that a does not change overly 
much when the interface geometry changes. Additionally, 
some fraction of the counterions will escape the corona, 
an effect we also ignore in our approximations. This 
will lead to region near the edge of the brush which is 
non-neutral. However, we expect that the size of such 
a region will be proportional to the electrostatic screen- 
ing length, which is small compared to the brush size. It 
seems doubtful that these approximations will change the 
overall qualitative conclusions of this paper, though they 
are likely to play a quantitative role in some aspects of 
our calculation. This will be discussed where appropriate 
in the following. 

Given the energetics described here we obtain mean- 
field results by minimizing the free energy per chain over 
unconstrained parameters. For the case of the equilib- 
rium calculations discussed below, we minimize all free 
parameters. The counterion free energy density does not 
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FIG. 3: The phase behavior for hairy micelles (/ = 0.3). The 
free energy per chain of cylinders and individual spheres is 
equal along the dashed line. The sphere phase is stable in 
the lower right portion of the phase space, that is, at high 
a and low Cg. Deep in the regime in which spheres have the 
lowest free energy, the cylinder state becomes unstable to the 
formation of undulated cylinders; this region is surrounded 
by the solid line. 



depend on the core size i?c, allowing us to minimize over 
Rc analytically, only on the ratio of brush height to core 
size, Tib- We retain three parameters to be minimized 
over numerically: the brush height /i, the undulation pa- 
rameter e, and the brush chain direction p. This gives us 



2^2 



po{^(eW8 + 3/S(H6,e,/3)/[^(l 
where / = 7V,/(7V, + TV,). 



/)]}' 



(13) 



II. RESULTS 
A. Metastability of Undulated Cylinders 

We computed mean-field phase diagrams for "crew- 
cut" (/ = 0.3), symmetric (/ = 0.5) and "hairy" 
(/ = 0.7) micelles. The free-energy minimization was 
performed for fixed charge fraction, a, and salt concen- 
tration, Cg. We assume for simplicity that po = and 
we consider copolymers composed of N = Nc-\- = 200 
monomers per chain, roughly equivalent to the experi- 
ments of Bendejacq et al 13]. We set the surface tension 
to ^o? = bhsT and find that hairy micelles (/ = 0.3) are 
computed to have core radii in the range of 5 to 10 nm 
(assuming a = 5 A) , which is in good agreement with the 
experimental observations of ref. 13]. 

The phase behavior for hairy micelles (/ = 0.3) as a 
function of charge fraction a and dimensionless salt con- 
centration Cga^ is shown in Fig. [HI The free energy per 
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FIG. 4: Plot of e and as a function of a for hairy micelles 
(/ = 0.3) at salt concentration Csa^ = 10 As a increases 
[3 is decreasing, indicating that the chains are tipping up ra- 
dially. 



chain of cylindrical and spherical micelles are equal along 
the dashed line. At high salt and low polymer charge 
cylindrical micelles are stable, while at low salt and high 
charge spherical micelles become stable. However, cylin- 
drical micelles that persist into the stability region of 
spheres become unstable to undulation. These undula- 
tions are driven by the frustration induced by the counte- 
rions, which prefer interfaces with positive Gaussian cur- 
vature in order to increase their free volume. However, 
this tendency toward the formation of spheres must be 
balanced by the significant increase in free ion concentra- 
tion in the regions of negative Gaussian curvature which 
is only partially mitigated by the migration of chains to- 
ward the positive curvature region (i.e. surface density 
of chains is reduced in the saddle regions; see Fig. [2 for a 
schematic description). Therefore we conclude that the 
undulated cylinder shape is metastable. 

In Figure we plot e and /3 as function of a for fixed 
salt concentration (c^a^ = 10~^). When the cylinders 
first become unstable at low charge fraction, the undu- 
lated cylinders have eccentricities on the order of e ^ 0.3 
and ~ 1 so that the chains are almost nearly directed 
normal to the interface. As the charge fraction is in- 
creased, the increased osmotic pressure within the brush 
increases its height (see Fig. 0|. This effect leads to 
an increase in the monomer and counterion concentra- 
tion in the negative curvature regions. The subsequent 
rise in osmotic pressure leads to a tendency of the chain 
trajectories to tip toward the radial direction. Once P 
decreases to zero and the brush chains extend radially 
from the cylinder axis, the counter ions no longer gain 
any entropy from the extra interfacial curvature of the 
unduloid. In fact, the core free energies, the core stretch- 
ing and interfacial energy, are minimal for the e = 
cylindrical state. And since the brush free energy is not 
lower for an unduloid state with /3 = than the e = 
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FIG. 7: The end of a cylinder is a piecewise continuous con- 
catenation of regions with constant interfacial curvature. A 
hemispherical cap is attached to a pair of unduloids which are 
attached to a cylinder. 



FIG. 5: The phase behavior for symmetric micelles (/ = 0.5). 
The free energy per chain of cylinders and individual spheres 
is equal along the dashed line. Here, the unduloid state is 
met ast able below and to the right of the dark line. 
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FIG. 6: The phase behavior for crew-cut micelles (/ = 0.7). 
The free energy per chain of cylinders and individual spheres 
is equal along the dashed line. Here, the unduloid state is 
met ast able below and to the right of the dark line. 



state, e decreases back to for large a. The boundary of 
the undulated cylinder region at high charge fraction is 
determined by the condition e = (see Fig.0J. 

The phase behavior of symmetric (/ = 0.5) and crew- 
cut (/ = 0.7) micelles is shown in Figs. and El re- 
spectively, and is very similar to the hairy case. As 
we progress toward the crew-cut limit, the CS transition 
moves toward larger charge fractions and lower salt con- 
centrations. This may be due to the core energy which 
prefers cylindrical micelles over spheres and is dominated 
by the interfactial term, {TJ, becomes more important as 
the size of the core block, Nc^ increases. As the core 



chain length increases, the undulated cylinder behavior 
no longer disappears at large charge fraction for a < 1. 
As in the case of hairy micelles, however, spheres become 
energetically favorable well before the UC region in both 
cases. 



B. End caps 

Given our model for the free energy of micelles con- 
structed from the unduloid surfaces, we build an ansatz 
for the shapes of the ends to micelles of cylindrical mi- 
celles from consider piecewise constant mean curvature 
surfaces. This is done by "gluing" two unduloids quarter- 
periods and a spherical hemisphere to the end of a cylin- 
der with radius Rc (see Fig. Q. One of the unduloid 
pieces changes radius from Rc to Rn^ the neck radius. 
The second piece increases the radius from Rn to the 
spherical cap radius, Rg. Since we assume that the mi- 
celles are effectively infinite, the end pieces can increase 
in volume by pulling chains from the infinite cylindrical 
portion, and thus, raise or lower the energy of that par- 
ticular chain relative to its energy in the cylinder. Thus, 
to compute the free energy of these ends we consider that 
the ends are in contact with an infinite source of poly- 
mer chains at fixed chemical potential, which is set by 
the free energy per chain of the infinite cylindrical state. 
The cylinder radius, Rc^ is set to the equilibrium value 
of an infinite cylinder, and we minimize the free energy 
of the end over all remaining unconstrained parameters. 

Many of the features of the characteristic end shapes 
can be explained by considering the free energy per chain 
in different geometries. In the region of cylinder stability 
(at low charge and high salt), chains in the endcap have 
a higher free energy than those in the bulk cylinder, and 
therefore preferentially migrate into the cylinder region. 
This has the effect of decreasing the radius of the spher- 
ical endcap, giving the appearance of the interface 
at a = 0.1 in Fig. |H1 In the spherical stability regime 
(at high charge and low salt), we find that the endcap, 
i?s, is generically larger than the cylinder radius, Rc. 
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FIG. 8: Typical end shapes for various charge fractions, in 
particular, for hairy micelles (/ = 0.5) at salt concentration, 
Cstt^ = 10~^. The cylinder radius Rc is set by the equilibrium 
cylinder radius. 



Since in this regime the free energy of chains within the 
spherical bulges is lower than for chains with the cylin- 
drical portion, chains preferentially migrate toward the 
spherical cap, increasing the radius of end cap, Rg. To 
best approximate the preferred spherical geometry, the 
neck radius i?^, pinches down to well below Rc. This 
procedure yields a spherical cap radius Rs very close to 
the ideal spherical micelle size. This characteristic "bul- 
bous" end shape has been widely observed in a number 
of polymeric micelle systems 0, , and is similar to 
the ends conjectured to appear in the context of bicon- 
tinuous phases of emulsions 25^. However, despite the 
preference for a spherical micelle configuration indicated 
by these bulbous end shapes, the free energy always has 
a minimum at finite neck size. Indeed, the high osmotic 
pressure induced by the large ion concentration in the 
saddle regions of the charged brush domain lead to siz- 
able energy barriers to "pinching off" {Rn 0). These 
barriers are discussed in detail below. 



C. Energy Barriers to CS Transition 

Since undulated cylinders are a metastable state, we 
consider here the free energy barriers governing the even- 
tual breakup of the micelles into spheres. Diblock copoly- 
mer kinetics is a subject of ongoing theoretical [2y, [23 
and experimental [2^ research. Since diblocks in sol- 
vent are the polymeric equivalent of amphiphilic surfac- 
tants, it is instructive to compare the kinetics of their 
respective micellar distributions. In surfactants, there 
is a substantial fraction of free molecules and the pri- 
mary kinetic mechanism for morphological changes is the 
Aniansson-Wall (AW) mechanism [sO, |^. Briefly, the 
exchange of single surfactant molecules (unimers) from 
micelle to bulk solution is fast enough to allow the di- 
rect formation of spheres and the subsequent depletion 



of molecules from cylindrical to spherical micelles. In 
the higher molecular weight copolymers, however, the 
concentration of free chains is significantly suppressed 
and, consequently, the AW unimer exchange mech- 
anism is drastically slowed 29]. Instead, Dormidontova 
has found that fission-fusion events of micelles likely play 
a significant role in the relaxation kinetics of copolymers 
solutions 27]. This suggests that the kinetic pathway 
connecting cylindrical and spherical micelle distributions 
may be dominated by states of fixed polymer number 
and intermediate interfacial geometry, such as undulated 
cylinders |9|. 

We turn our attention to the calculation of the free en- 
ergy barriers for the fission of a cylindrical (or undulated) 
micelle end into spherical micelles. Our end shape ansatz 
allows us to compute the energetic barrier for budding a 
sphere from the end of a cylindrical micelle by continu- 
ously tuning the neck size, Rn. A schematic free energy 
plot for this pathway is depicted in Fig.[Ul We consider a 
pathway which is divided into three stages: (1) contrac- 
tion of the neck radius, (2) scission of the neck, and (3) 
bud relaxation. While, in principle, we could imagine a 
pathway in which Rn contracts continuously to zero (or 
at least some molecular size), we find that allowing for 
neck scission lowers the over barrier considerably. There- 
fore, once the neck size has decreased below some critical 
size, the bud neck breaks, thereby exposing an area of the 
core 27ri?^ to solution. At the point of neck scission, we 
estimate the cost of exposing core monomers to solution 
as 7, the interfacial surface tension, times the exposed 
surface area. To compute the minimum energy pathway 
within our ansatz we minimize the free energy barrier, 
AFi)ar^ over this critical neck radius. 

Our calculation for the barrier height for end fission of 
spherical micelle from a cylinder end is plotted in figure 
El for hairy micelles (/ = 0.3) and for low salt (c^a^ = 
10~^). We find that AF^ar decreases for increased a. 
When a is small, the barrier heights can be quite large 
(thousands of ksT). For larger a, however, the barrier 
heights fall precipitously to tens or hundreds of kuT. In 
comparison to the case neutral copolymer melts |33 , the 
free energy barriers computed here are quite large, up to 
1000 times larger. 

A simple consideration of the free energy scale set by 
the brush free energy shows that the excluded volume 
interaction induced by the brush electrostatics should in- 
deed lead to free energy barriers of this magnitude. While 
excluded volume effects are screened in molten polymer 
brushes, the CS transition is frustrated by the increased 
osmotic pressure in the neck region of the end bud. Typ- 
ically, such an end bud will contain 10^ — 10^ polymer 
chains, and each of those chains entrains aNi) counter 
ions within the dissolved brush. Given that each of those 
counterions represents roughly ksT worth of free energy 
due to the entropy lost in confinement, it is then not sur- 
prising that such a topological transition could cost thou- 
sands of ksT. Nonetheless, we see from Figure 11171 that 
changing the effective charge of the polyelectrolyte block 
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AF 




FIG. 9: A skematic plot of the free energy along the assumed 
kinetic pathway for cylindrical ends pinching off to for spher- 
ical micelles. The horizontal axis is some generalized reaction 
coordinate which carries the micelle end through the CS tran- 
sition along the lowest energy pathway. The transition occurs 
in three stages: neck contraction below the equilibrium size, 
Rn, neck scission, and bud and and end relaxation. Presum- 
ably, the scission must expose the insoluble blocks to solvent, 
and therefore it is visualized as a discontinuous jump in the 
free energy of the end. 



^Ftar 




a 



FIG. 10: A plot of the CS barrier height as a function of a 
for hairy micelles (/ = 0.3) at low salt, Csa^ — 10~^. The 
inset displays the small barrier height region at large charge 
fraction in rescaled coordinates for convenience. 



can decrease the barrier at least two orders of magnitude. 
The qualitative explanation of this effect is simple. As a 
increases, the equilibrium size of the micelles decreases, 
indicating that fewer chains comprise the micelle ends, 
leading to an overall decrease in AF^ar- 

The barrier heights for hairy micelles are plotted on a 
log-log scale in figure^Jfor a number of different salt con- 
centrations. We find two characteristic regimes of power 
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FIG. 11: The CS barrier height for hairy micelles (/ = 0.3) 
as a function of a on a log-log scale for a number of different 
salt concentrations. 



law behavior, one at low charging and one at high charg- 
ing. At high salt concentration and large a, we find that 
the barrier height decreases as . In this regime, the 
decrease in neck radius Rn is small and the free energy 
barrier is dominated by direct scission. Here, the depen- 
dence on a comes directly from the equlilibrium scaling 
of the cylinder radius with a at high salt in the limit 
H}) i?c, Rc ^ OL~^^'^ Largely, we expect the equi- 
librium neck radius to scale with Rc^ Therefore, direct 
scission predicts a free energy barrier, AFi^ar ^ T^n? 
which yields the observed power law. 



III. DISCUSSION 

We have developed a free energy expression for di- 
block copolymer cylinders with a hydrophobic block and 
a charged corona, and used this to understand the stabil- 
ity of various axisymmetric structures. We use a family 
of constant mean curvature surfaces called unduloids to 
parametrize possible interfaces from cylinder to a neck- 
lace of spheres. This ansatz allows us to probe the 
full, non-linear dependence of the free energy on the mi- 
celle surface geometry. Our results indicate that there 
is a metastable undulated cylinder phase, as observed in 
experiments J^J which occurs generically in parameter 
ranges where the free energy per chain for spherical mi- 
celles is lower still. We have studied the shape of the 
spherical endcaps and found, in particular, that in the 
regime in which cylinders are stable, the ends show a 
distinct decrease in radius below the equilibrium value. 
When spherical micelles are preferred in equilibrium, the 
spherical cap radius is larger than the cylinder radius - 
and is close to the preferred spherical micelle size. We 
use these results to compute dependence the free energy 
barriers for budding spherical micelles from cylindrical 
micelles on polymer charge, a, and solution salt concen- 
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tration, c^. 

Our calculation of the free energy is variational and 
therefore our values for the free energy barrier to bud- 
ding, AF^cir, necessarily provide upper bounds. However, 
certain conclusions from this study are robust. Indeed, 
these results along with the fact that the AW step-wise 
kinetic mechanism is suppressed in polymeric systems 
[27] strongly support the suggestion that distributions 
of polymer micelles in selective solvent "non-ergodic," 
with micelles essentially trapped within a fixed, quenched 
topology 9]. In contrast with neutral copolymer solu- 
tions, however, we find here a strong, tunable dependence 
of barrier height on polymer charge and salt concentra- 
tion. Kramers kinetic theory predicts that the CS tran- 
sition rate, kcs should depend on the the free energy as 
kcs cx e-^^^--/^^^ At smah a, the CS pathway 

investigated here is kinetically arrested by large barrier 
heights. Supposing the AW unimer exchange mechanism 
is also sufficiently slow, these large barrier heights suggest 
that the metastable undulated cylinder regime may be 
observed. Such undulated cylinders have been observed 
experimentally in a number of amphiphilic polymeric sys- 
tems, both charged and uncharged j9|, [lO|, IjJ • At large a 
and low Cs , AFi^ar may fall to sufficiently low values that 
the transition to spherical micelles may proceed on ex- 
perimental time scales. This suggests that the observed 
"phase boundary" between cylindrical and spherical mi- 
celles in ref. 13], may reffect that fact that these kinetic 
barriers depend quite strongly on a and Cg. Given that 
AFi)ar varies by two orders of magnitude, we expect an 
extremely sharp dependence of the kcs on a and Cg over 
the same parameter range. 

The geometry of self-assembled structures is often cast 
in terms of constant mean curvature (CMC) surfaces |25|, 
I34] . More generally, the energy of any interface can be 
expanded in powers of the curvature, resulting in the 
Helfrich energy. 



duloid family of surfaces capped at both ends. The en- 
ergy reduces to 



E = -fAt 



dA {H-Ho 



dA K, (14) 



where Atot is the total micelle interfacial area, H = 
(l/2)(l/ri + 1/^2) is the mean curvature, K = l/(rir2) 
is the Gaussian curvature, Hq is the spontaneous curva- 
ture, which comes from the molecular asymmetry of the 
chains and the corona-solvent interaction. The parame- 
ter 7 is the surface tension. Finally, K: and R are effective 
bending rigidities for the mean and Gaussian curvatures 
respectively, which must be derived from a microscopic 
model. For = ^ = 0, this energy describes the shape of 
soap bubbles, and CMC surfaces minimize the energy in 
this limit. In the context of diblock copolymers 1^ 
or ffuid membrane shapes the additional curvature 
terms modify the shape and stability of various struc- 
tures. Higher order curvature terms are typically ne- 
glected as being subleading, but we will show equation 
(|T^ fails to describe our results, even qualitatively, with- 
out them. 

We restrict ourselves to interfaces belonging to the un- 



E 



7- 



Atot{e) - ^TTi^, (15) 



where we note that a cylinder with two ends has 
J dA K = — 47r from the Gauss-Bonnet theorem. If we 
choose the length of the cylinder to conserve the vol- 
ume, the total area is given by Atot{^) = [A{e)/V{e)]Vtoti 
where A{e) and V(e) are the reduced area and volume of 
a half undulation wavelength (see appendix) and Vtot is 
the volume of the cylinder. The only dependence on e 
comes from the volume and area of the interface (recall 
that H does not change with e) which, at fixed Rc^ is 
minimum for e = 1 no matter the value of spontaneous 
curvature. This refiects the fact that spheres minimize 
the area for a fixed volume, and indeed, this is a man- 
ifestation of the well-known Rayleigh instability of fiuid 
interfaces 38]. The spontaneous curvature only sets the 
overall size scale of the micelles, not their topology. 

We conclude that higher order curvature terms to 
equation ([T^ are required to understand the local sta- 
bility of undulated cylinders in terms of universal geo- 
metrical measures of the block interface (2^|3^. This is, 
perhaps, not surprising since higher order corrections in- 
troduce length scales on the order of Rc for axisymmetric 
shapes In the case of vesicles, where equation (|T^ 
is known to apply, these corrections are controlled by the 
membrane thickness [s^, 11^ • We leave the development 
of a geometrical shape equation for diblock copolymer 
micelles to future work. 
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APPENDIX A: UNDULOID GEOMETRY 

In this Appendix we present details concerning the ge- 
ometry of unduloid surfaces which will be useful in our 
free energy calculation. Given the parametric represen- 
tation of X(ii, (/)) by eqs. ^ and Q we have the differ- 
entials. 



dr = 



Rc{e cos u)du 



2e sinu 



(Al) 
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and 



dz 



i?c(l + e sin u)du 
Vl + e2 + 2e 



(A2) 



The tangent vectors for the two principle directions of 
the surface are given by, 



(iX^ = r{u 



(A3) 



and 
(iX, 



Rcdu 



: [{1 -\-e sin u)z-\-e cos ur] . (A4) 



Vl + e2 + 2e 

From these we compute the differential element of area, 



2esinudud(j) , (A5) 

and the surface normal, 

dX(f) X dXu [- ecosuz {1 esinu)f] 



\dXs X dX^I 



Vl + e2 + 2ee 



(A6) 

Since the surface is periodic we concern ourselves only 
with a half period of the undulation, u G [— 7r/2, 7r/3]. 
The volume per half period can be computed, 

r^^ dz 

V(e) = TT / du^r\u) 
' du 



7r/2 



/ duy 

J-7r/2 



1 + + 2e sin u{l + e sin u) . 

(A7) 



The dimensionless volume, V(e), is then defined by 
V{e) = 7r^i?^V(e). Finally, we surface area of a half pe- 
riod is computed by, 

^ 

A{e) = 2iiRl / du^l^e^^ 2e sin u , (A8) 

and the dimensionless area, w4(e), is defined by A{e) = 
2^^RlA{e). 



APPENDIX B: MICELLE CORE ENERGETICS 

In this section we compute the core stretching (0) . 
As described in Sec. |D we divide the core region into 
infinitestimal wedges, which contain polymer chains ex- 
tending in the same direction. We assume that the core 
chains (and wedges) extend along the surface normal. As 
shown in Figure IT^ the vector, T>c{u) = r{u)r — j^r{u)z^ 
describes the extension of the wedge. From the normal 
stretching condition, we have the mapping from a point 
on the surface to a point along the z axis, located a zo{z) 
(see Figure [T^ . This mapping is given by. 



dv 

zoiu) = z{u) + ■ 



(Bl) 




zo{u) Zs{u) 

FIG. 12: A schematic representation of the normal trajectory 
ansatz. 



In general, this mapping must be consistent with the con- 
straint the chain paths (and wedges) do not cross. This 
means that the mapping must be one-to-one. This con- 
straint is satisfied for ^ > 0. This yields an inequality 
which must be satisfied by surfaces we consider. 



d^r /dr\'^ 
d^^'^Uz) 



< 1 



(B2) 



It is not difficult to show that all unduloid surfaces satisfy 
this inequality. 

To compute the stretching integral defined in eq. ©, 
we must compute the volume distribution within the 
wedge, which we do by considering the dimensions of the 
wedge in the d'Ku and dX^directions as a function of dis- 
tance to the z axis, h^. We will denote these directions 
by d£u{h,u) and d£(f){h,u), respectively (see Figure [T3|) . 
These are given by, 

d£u{ho,u) =dzo{n-r)^ J^^ {Rcdu-dzo{n-r)) , (B3) 



and 



di^{hc,u) 



ho 



Dc{u) 



r{u)d(j) 



(B4) 



where Dc{u) = \Y>c{u)\. From this we compute the 
area distribution in the core wedges by dAc{ho^u) = 
diu{ho^u)di(/y{ho^u). With this the volume per core 
wedge can be computed. 



dvc{u) = / dhodAc{ho^u) 
Jo 



(B5) 



and the surface density, cr(ix), of chains reaching the in- 
terface at. 



cr{u) 



dVc{u)po 
NcdA{u) ' 



(B6) 



We plot a{u) for various unduloid shapes in Figure 
For cylinders (e = 0) and spheres (e = 1) the surface 
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Hb{d-bb-Tb)+\dXu\{n-t)i,) 



dzo{n-r) 

FIG. 13: A differential wedge extending along the mean chain 
trajectory. 




cr(M,0.15) 




FIG. 14: Plots of the surface chain density for various undu- 
loid shapes. The density is plotted in units of poRc/Nc. 



density is constant, while for intermediate shapes a{u) 
drops in the saddle regions. 

Finally we compute the stretching free energy for a half 
period of the core region, 



QAT2~^ / # / du dho 
8N^a^ Jo J-7r/2 Jo 

x{Dc{u) - hofdAc{ho,u) 



■S{e) , 



(B7) 



where the dimensionless stretching is defined on the last 
line. Dividing this with the number of chains in a half 
period of core is 7r^i?^V(e)po/^c, we obtain eq. ©. 



APPENDIX C: OUTER BRUSH ENERGETICS 

In Figure we show the geometry of the chain tra- 
jectories in the outer polyelectrolyte brush. The vec- 
tor D5(i^) = Hi){u)'Db{e^ /3) gives the extension of brush 
chains from the interface surface at X(i^, cf)) to the sur- 
face described by the tips of this outer brush at X(ix, cj)) + 



FIG. 15: Geometry of coronal polymer brush. 



^^{u). As described in Sec. ^make the assumption that 
all chains are uniformly stretched, with their free ends 
constrained to tips of the brush. Given this assumption, 
it is straightforward to compute the free energy cost of 
extending the brush within a half period of undulation. 



7r/2 

du I duH^{u)a{u) 



dA 
du 



(CI) 



Dividing by the number of chains in the core, we obtain 
eq. (Uni). 

Next we compute the brush monomer concentration, 
Cm(r, which depends on the brush geometry. There- 
fore, we need to determine the volume distribution as 
we did for the core brush. In particular, we divide the 
outer brush into wedges extending along the mean poly- 
mer trajectory and we compute the area distribution 
within these wedges as we did in the case of the core (see 
Fig. [T5|) . Here, we define ^^^(ii, as we did before, 
and diu{h^hi)) is the dimension of the wedge associated 
with the direction perpendicular to D5 which we denote 
T5 = (D5 •r)z— (D5 • z)r. These can readily be computed 
from the geometry. 



dlu{K h) = \dXu\{h • D5) + hidBt ■ %) , (C2) 



and 



di^{u, h) = r{u)d(^ + /i5(D5 • r) . (C3) 



We can compute the area available to polymers in 
the wedge at a distance, from the interface by 
dAi){u^ h})) = d£u{u, hi))d£(f){u, hi^). If a given wedge leaves 
the surface from X(i^,(/)), then it contains a{u)dA{u) 
polymer chains. Given that we consider the corona chains 
to be evenly stretched throughout the brush, a chain 
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leaves, dhi^Ni^/Hiy monomers between and + We use this relation in eq. (|T^ to compute the concen- 

in a differential volume, dAi){u^hb)dhb. Therefore, the trations of free ions in the brush, 
position dependent monomer concentration is, 

Nba{u)dA{u) 
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